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The most energetic particles ever detected exceed 10 20 eV in energy. Their existence represents at the same 
time a great challenge for particle physics and astrophysics, and a great promise of providing us for a probe of 
the validity of the laws of Nature in extreme conditions. We review here the most recent data and the future 
perspectives for detection of cosmic rays at ultra-high energies, and discuss possible ways of using these data to 
test the possibility that new Physics and/or new Astrophsyics may be awaiting around the corner. 

"When you carry out an experiment there are two possible outcomes: 
either you confirm the theoretical expectation, and in this case you made 
a measurement, or you don't, and in this case you made a discovery" . 

E. Fermi 



1. Introduction 

The cosmic ray spectrum has been now mea- 
sured over a large range of energies, that extends 
in its upper part to more than 10 20 eV, the so- 
called ultra high energy cosmic rays (UHECRs). 
The quest for the origin of these high energy 
particles over the whole range is still open and 
represents one of the big challenges for the fu- 
ture. In the highest energy end of the cosmic 
ray spectrum, several issues make the challenge 
even harder: a) it is hard to envision possible 
acceleration sites where particles with energy in 
excess of 10 20 eV may be accelerated; b) even 
if some class of sources could indeed accelerate 
particles to the highest energies, a homogeneous 
spatial distribution of these sources would leave 
an imprint in the cosmic ray spetrum, known as 
Greisen-Zatsepin-Kuzmin (GZK) cutoff, due to 
photopion production on the photons of the cos- 
mic microwave background jjj. On the basis of 
the first issue, we would expect an end of the cos- 
mic ray spectrum to occur at some high energy. 
On the basis of the second issue, we would expect 
a strong flux suppression (not really a cutoff) at 
about 5 x 10 19 eV. 

Several experiments have been operating to de- 
tect the flux of UHECRs, starting with Volcano 



Ranch and continuing with Haverah Park || 
and Yakutsk [ 4 1 to the more recent experiments 
like AGASA |-|, Fly's Eye |L|n| and HiRes 

& 

At present there is no clear indication that the 
cosmic ray spectrum comes to an end due to ei- 
ther one of these two reasons. More statistics of 
events is however required to achieve a solid con- 
clusion in this respect. 

The physics involved in the explanation of the 
origin and propagation of UHECRs needs often to 
be pushed to its extremes to accomodate observa- 
tions. This transforms a problem into a precious 
tool to probe a territory which is still uncharted, 
a New Physics which several hints tell us should 
exist, but is still hidden somewhere. In the fol- 
lowing we describe a few directions in which this 
investigation may lead. 

In §[2| we briefly summarize the status of current 
observations; in §||wc restate the basic issues that 
are known as GZK problem; in §|] we illustrate 
two examples of new physics that can be investi- 
gated with the help of UHECRs, na mel y Physics 
close to the grand unification scale (§ |4.l|) a nd pos- 
sible violations of Lorentz invariance (^4^2). In §|5| 
we give our conclusions. 
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2. Observations 

The cosmic ray spectrum is measured from 
fractions of GeV to a (current) maximum energy 
of 3 x 10 20 eV. The spectrum above a few GeV 
and up to ~ 10 15 eV (the knee) is measured to be 
a power law with slope ~ 2.7, while at higher 
energies and up to ~ 10 19 eV (the ankle) the 
spectrum has a steeper slope, of ~ 3.1. At en- 
ergy larger than 10 19 eV a flattening seems to be 
present. 

The statistics of events is changing continu- 
ously: an analysis of the "all experiments" statis- 
tics was carried out in [H and found 92 events 
above 4 x 10 19 eV. 47 events were detected by 
the AGASA experiment. A more recent analy- 
sis H of the AGASA data, carried out expanding 
the acceptance angle to ~ 60°, has increased the 
number of events in this energy region to 59. 

In (|] the directions of arrival of the AGASA 
events (with zenith angle smaller than 45°) above 
4 x 10 19 eV were studied in detail: no apprecia- 
ble departure from isotropy was found, with the 
exception of a few small scale anisotropies in the 
form of doublets and triplets of events within an 
angular scale comparable with the angular reso- 
lution of the experiments (~ 2.5° for AGASA). 
This analysis was repeated in jl3| for the whole 
sample of events above 4 x 10 19 eV, and a total 
of 12 doublets and 3 triplets were found within 
~ 3° angular scales. The attempt to associate 
these multiplets with different types of local as- 
trophysical sources possibly clustered in the local 
supercluster did not provide evidence in that di- 
rection Q. 

Recently, the AGASA collaboration reported 
on the study of the small scale anisotropies in the 
extended sample of events with zenith angle < 
60°: 5 doublets (chance probability ~ 0.1%) and 
1 triplet (chance probability ~ 1%) were found. 

The information available on the composition 
of cosmic rays at the highest energies is quite 
poor. A study of the shower development was 
possible only for the Fly's Eye event jll] and dis- 
favors a photon primary []l5|| . A reliable analy- 
sis of the composition is however possible only 
on statistical basis, because of the large fluctu- 
ations in the shower development at fixed type 



of primary particle. The Fly's Eye collaboration 
reports evidence for a predominantly heavy com- 
position at 3 x 10 17 eV, with a smooth transition 
to light composition at ~ 10 19 eV. This trend was 
later not confirmed by AGASA |j 16 1. Recently 
in ref. JTjJ the data of the Haverah Park experi- 
ment on highly inclined events were re-analyzed: 
this new analysis results in no more than 30% of 
the events with energy above 10 19 eV being con- 
sistent with photons or iron (at 95% confidence 
level) and no more than 55% of events being pho- 
tons above 4 x 10 19 eV. 

Recently new data have been presented by the 
HiRes experiment |TJ] . These new results seem to 
be in agreement with the presence of a GZK fea- 
ture in the cosmic ray spectrum, and are therefore 
in disagreement with the AGASA data. In order 
to quantify the discrepancy we ran some simu- 
lations jlTJ of the statistics of events expected 
with the AGASA and HiRes exposures. In table 
I we report our results: the first column gives the 
energy threshold used to measure the number of 
events quoted in the other columns. The injection 
spectrum is taken as E~ 2A , because it provides 
the best fit to the data points of AGASA at lower 
energies. The evolution in the luminosity of the 



sources is taken as L cx (1 



\m+3 



where the 3 



accounts for the redshift (to = corresponds to 
no intrinsic evolution). For AGASA we also con- 
sidered the possibility that there is a systematic 
error of 10% in the energy determination (column 
indicated as m = + 10% sys.). From the num- 
bers in the table we deduce that the AGASA data 
do not agree with the GZK prediction at the level 
of 2.4cr (3.6(7) for m = (to = 4). On the other 
hand HiRes data seem to agree with the predicted 
flux. 

Clearly we need more data in order to finally 
settle the issue. 

3. The physical meaning of the GZK fea- 
ture 

The puzzle of UHECRs can be summarized in 
the following points: 



The production problem: 



Acceleration 



mechanisms have to be pushed to their 
extremes in order to allow the production 
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Table 1 

Simulated events for the AGASA and HiRes exposures. The injection spectrum is taken as E~ 2A . The 
source evolution is taken proportional to (1 + z) 3+m . 



of particles with energy in excess of 10 20 
eV. 

• The large scale isotropy: observations show 
a remarkable large scale isotropy of the ar- 
rival directions of UHECRs, with no cor- 
relation with local structures (e.g. galactic 
disk, local supercluster, local group). 

• The small scale anisotropy: the small (de- 
gree) scale anisotropics, if confirmed by fur- 
ther upcoming experiments, would repre- 
sent an extremely strong constraint on the 
type of sources of UHECRs and on mag- 
netic fields in the propagation volume. At 
present no association has been found be- 
tween the arrival directions of UHECRs and 
nearby bright sources. 

• The GZK feature: the GZK cutoff is mainly 
a geometrical effect: the number of sources 
within a distance that equals the pathlength 
for photopion production is far less than the 
sources that contribute lower energy parti- 
cles, having much larger pathlength (com- 
parable with the size of the universe). The 
crucial point is that the cutoff is present 
even if plausible nearby UHECR engines are 
identified. 

The solution of the problem of UHECRs in its 
several aspects will come when many pieces will 
fall into place. A careful measurement of the 
spectrum of UHECRs and of their chemical com- 
position will play a crucial role for the identifica- 
tion of the sources. In particular, the composition 
may be a smoking gun in favor or against whole 
classes of models. 

In this section we consider in some more detail 
the issue of the GZK feature and the implications 



of its possible absence, suggested by the AGASA 
data. 

It is often believed that the identification of 
one or a class of nearby UHECR sources would 
explain the observations and in particular the ab- 
sence of the GZK feature. This is not necessarily 
true. The (inverse of the) lifetime of a proton 
with energy E is plotted in fig. 1 (left panel) 
together with the derivative with respect to en- 
ergy of the rate of energy losses b(E) (right panel) 
[the figure has been taken from ref. pcfl l. The 
flux per unit solid angle at energy E in some 
direction is proportional to n \(E)&(E), where 
no is the density of sources (assumed constant), 
X(E) = c/((l/E)dE/dt) and $(£) is the source 
spectrum. This rough estimate suggests that the 
ratio of detected fluxes (multiplied as usual by 
E 3 ), at energies E\ and E 2 is 

_ EfF(Ex) XjE^jE^E 3 _ AQgi) f E ± \ 
E$F{E 2 ) ~ \(E 2 )<5>{E 2 )E 3 \{E 2 ) \E 2 ) 

where in the last term we assumed that the source 
spectrum is a power law §(E) ~ _E -7 . If for in- 
stance one takes E\ = 10 19 eV (below Eqzk) and 
E 2 = 3 x 10 20 eV (above Eqzk), from fig. 1 one 
obtains that 1Z ~ 80 for 7 = 3 and 1Z ~ 10 for 
7 = 2.4. The ratio 1Z gives a rough estimate of 
the suppression factor at the GZK cutoff and its 
dependence on the spectrum of the source. For 
flat spectra (7 < 2) the cutoff is less significant, 
but it is more difficult to fit the low energy data 
|||] (at E ~ 10 19 eV). Steeper spectra make the 
GZK cutoff more evident, although they allow an 
easier fit of the low energy data. The simple ar- 
gument illustrated above can also be interpreted 
in an alternative way: if there is a local overden- 
sity of sources by a factor ~ 1Z, the GZK cutoff 
is attenuated with respect to the case of homoge- 
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Figure 1. From jf^j/. Left panel) (l/E)dE/dt for 
a proton in JjZ|/ ( curve 1 ), in Jjl|/ ( curve 2) and in 
Jjj|/ (black squares). The curve 3 is the contribu- 
tion of the red shift. Right panel) The derivative 
db(E)/dE, with b(E) =dE/dt at z = 0. 



neous distribution of the sources. The question of 
whether we are located in such a large overden- 
sity of sources was addressed, together with the 
propagation of UHECRs, in Assuming that 
the density of the (unknown) sources follows the 
density of galaxies in large scale structure surveys 
like PSCz |Q and Cfa2 (25), the authors estimate 
the local overdensity on scales of several Mpc to 
be of order ~ 2, too small to compensate for the 
energy losses of particles with energy above the 
threshold for photopion production. 

As discussed in the current observational 
situation is not crystal clear in pointing toward 
the absence or the presence of the GZK feature in 
the data. Recalling the sentence of Enrico Fermi, 
quoted in the beginning of this paper, we are in 
a ironic situation in which we have a 2.5a discov- 
ery from AG AS A, and an equally significant mea- 
surement from the HiRes collaboration. What is 
absolutely clear is the need for higher statistics 
experiments, such as Auger and EUSO, that for- 
tunately are on the way. 

The bottom line of this section can be summa- 
rized in the following few points: 



1) the GZK cutoff is not avoided by finding 
sources of UHECRs that lie within the pathlcngth 
of photopion production, unless these sources are 
located only or predominantly nearby and are less 
abundant at large distances. 

2) To establish without any doubt the statis- 
tical significance of the GZK feature in the cos- 
mic ray data we need a statistics of events with 
energy > 10 20 eV larger than the present one. 
The enhanced statistics will also make easier to 
study the chemical composition and to identify 
the sources of UHECRs. 

4. Probing new physics 

UHECRs provide a laboratory where we may 
investigate the extremes of the Physics we know 
and possibly check whether there is something 
more than that. In this paper, we discuss only 
two examples, one related to the production of 
UHECRs (Top-down models) and the other, pos- 
sible minuscule violations of Lorentz invariance, 
that might affect the propagation of UHECRs 
on cosmological distances. Many more examples 
may be found but will not beb discussed here. A 
come comprehensive review can be found in |5q] . 

4.1. Top-Down (TD) models 

The old question of how the UHECRs are pro- 
duced is as fashionable now as it was 30 years ago. 
The severe requirements needed for ordinary ac- 
celeration processes to achieve ultra-high energies 
have fueled much interest in production mecha- 
nisms that work from the top, meaning that par- 
ticles are not accelerated but rather injected as a 
result of the decay of very massive unstable relics 
of the big bang. This may occur either as a result 
of decay processes of topological defects or by de- 
cays/annihilations of supermassive relic particles. 

Topological defects are naturally formed at 
phase transitions and their existence has been 
proven by direct observations in several experi- 
ments on liquid crystals and ferromagnetic ma- 
terials. Similar symmetry breakings at particle 
physics level are responsible for the formation of 
cosmic topological defects [for a review see [p8| ]. 

The fact that topological defects can gener- 
ate UHECRs was first proposed in the pioneering 
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work in Ref. [[29|. The general idea is that the 
stability of the defect can be locally broken by 
different types of processes and result in the false 
vacuum, trapped within the defect, to fall into the 
real vacuum (outside universe), so that the gauge 
bosons of the field trapped in the defect acquire 
a mass mx- At this point, the very massive and 
unstable particles rapidly decay producing high 
energy particles. 

Several topological defects have been studied 
in the literature: ordinary strings [^0|, supercon- 
ducting string s p9f , bound states of magnetic 
monopoles [[H 32 1, networks of monopoles and 
strings Q , necklaces |34) and vortons |35| . 

The basic idea, common to all TD models, 
is that the decay of a supermassive particle re- 
sults in the production of a quark-antiquark pair 
that hadronizes into mesons and protons. At the 
source, the composition of the produced parti- 
cles is dominated by gamma rays and neutrinos, 
while only about 5% of the energy goes into pro- 
tons. After propagation over cosmological dis- 
tances, the relative abundance of gamma rays and 
protons changes, as illustrated in fig. 2. Gamma 
rays become then dominant only at extremely 
high energies, although an appreciable fraction 
of the composition is still made of gamma rays. 
Much discussion exists on which topological de- 
fects may generate the observed fluxes, as sum- 
marized in [ ^6| , p7| . 

An interesting alternative to topological defects 
is represented by super-heavy (SH) relics of the 
big bang. SH particles with very long lifetime can 
be produced in the early universe and generate 
UHECRs at present ^f^Q. In the following, 
we will call these particles X-particles. 

X-particles can be produced in the early uni- 
verse through different mechanisms. The simplest 
of them is the gravitational production: particles 
are produced naturally in a time variable gravi- 
tational field or indeed in a generic time variable 
classical field. In the gravitational case no ad- 
ditional coupling is required (all particles inter- 
act gravitationally) . If the time variable field is 
the inflaton field 0, a direct coupling of the X- 
particles to <j> is needed. 

The gravitational production of particles was 
first proposed in |)43|. It does not require any 
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Figure 2. Spectrum of protons and gamma rays 
from topological defects ^o\. The two sets of 
curves for gamma rays are obtained for two dif- 
ferent estimates of the radio background. Differ- 
ent curves refer to different values of the mass of 
X-particles. 



additional assumption neither on the X-particles 
nor on cosmology. In particular inflation is not 
required a priori, and indeed it reduces the effect. 
It can be shown that at time t, gravitational pro- 
duction can only generate X-particles with mass 
mx < H(t) < m<f,, where H{t) is the Hubble con- 
stant and TO^ is the inflaton mass. It was demon- 
strated in (4^^l[l that the fraction of the critical 
mass contributed by X-particles with mx ~ 10 13 
GeV produced gravitationally is fix ~ 1, with 
no additional assumption! In other words, cold 
dark matter can naturally be explained in terms 
of X-particles in this range of masses. 

As mentioned in the beginning of this section, 
in order for X-particles to be useful dark matter 
candidates and generate UHECRs they need to 
be long lived. The gravitational coupling by it- 
self induces a lifetime much shorter than the age 
of the universe for the range of masses which we 
are interested in. Therefore, in order to have long 
lifetimes, additional symmetries must be postu- 
lated: for instance discrete gauge symmetries can 
protect X-particles from decay, while being very 
weakly broken, perhaps by instanton effects |$9|. 
These effects can allow decay times larger than 
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Figure 3. Spectra of gamma rays from the decay of 
SH relics. The parameters are chosen as indicated 
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the age of the universe, as shown in |44|| . 

The slow decay of X-particles produces UHE- 
CRs. The interesting feature of this model is 
that X-particles cluster in the galactic halo, as 
cold dark matter |||]. Hence UHECRs are ex- 
pected to be produced locally, with no absorp- 
tion. As a consequence, the observed spectra 
are nearly identical to the emission spectra, and 
therefore gamma rays dominate. The very flat 
spectra and the gamma ray composition are two 
of the signatures. The calculations of the ex- 
pected fluxes have been performed in p(i| . [l5| , [l6[ . 
In figure 3 we report a typical prediction of the 
spectra for this model, as derived in |3£|. The 
strongest signature of the model is the anisotropy 
due to the asymmetric position of the sun in the 
Galaxy [||]3||4|,|5(J . These papers agree that the 
present data are consistent with the anisotropy 
expected in the model of X-particles in the halo, 
due to the large uncertainties in the measurement 
of anisotropy. In fact, as discussed in §||, obser- 
vations at present do not show any appreciable 
deviation from isotropy, with the exception of a 
few degree scale anisotropics showing up in the 
form of doublets and triplets of events within an 
angular scale comparable to the resolution of the 
experiments. 

In TD models the presence of these multiplets 



of events is not well accomodated because of the 
homogeneous distribution of the topological de- 
fects or of the SH particles. However, it was 
shown in that in the latter the presence of 
the multiplets may actually be accomodated. 

As an alternative to the long lifetimes required 
by the decay of SH relics, a model based on their 
annihilation has been proposed |58), where how- 
ever large annihilation cross-sections are required. 

4.2. Lorentz invariance: is it really a good 
symmetry? 

UHECRs are the particles with the largest 
Lorentz factor that we have been able to mea- 
sure. They are then the most powerful probe to 
test the validity of Lorentz invariance (LI), one of 
the fundamental ideas of known Physics. LI viola- 
tions affect the propagation of UHECRs mainly 
because they may change the kinematic thresh- 
olds for photopion production, although other ef- 
fects have also been investigated. The pioneer- 
ing work in this direction was presented in |5q| . 
A very general parametrization of the violation 
of LI was discussed in |54j, where it was shown 
that the threshold may even move to infinite en- 
ergy of the ultra high energy proton, meaning 
that the reaction becomes kinematically forbid- 
den. The possible violation of Lorentz invariance 
at the Planck scale has been proposed by several 
authors as a possible explanation of the detection 
of events with energy above 10 20 eV. However, 
it is clear from the discussion in that the cur- 
rent number of such events is not sufficient to jus- 
tify such a proposal. Nevertheless, the shape of 
the spectrum of UHECRs can provide interesting- 
constraints on the strength of the LI violation, if 
the GZK feature is detected by upcoming exper- 
iments. 

Recently, an interesting extension of the previ- 
ous ideas on violation of LI was discussed in (5^] . 
In the new approach, the violation is induced by 
the quantum fluctuations of the space-time foam: 
if the Planck length is the smallest length which 
can be measured, then this scale also represents a 
lower limit on the error in a length measurement, 
and in particular of the Dc Broglie length of a 
particle. The induced uncertainty in the momen- 
tum of the particle is 8p ~ p 2 /Mp, so that energy 
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and momentum can be written as 



E = E 



E 2 



P 2 



— — p = p + Pirr- 



M 



Mr 



(2) 



where a and (3 are randomly distributed numbers 
with unit variance, and E and p are the average 
energy and momentum. Moreover the fluctua- 
tions in the metric can be written in the following 
form: 



Pn9 Pu 



2 , P 
P +7 Mp 



(3) 



where 7 is also a random number with normal dis- 
tribution and unit variance. When applied to the 
process of photopion production, the effect of the 
metric fluctuations is to move the threshold for 
photopion production around the classical value. 
However, as demonstrated in |56| ], when the reac- 
tion remains kinematically allowed, the threshold 
is moved to ~ 10 15 eV, so that on average the ef- 
fect of violation of LI due to a fluctuating metric 
is mostly of making low ener gy p articles vulnera- 
ble to photopion production |5g] 

5. Conclusions 

The search for the end of the cosmic ray spec- 
trum, started a few decades back in time, is still 
ongoing and still not successful. This challenge 
led us to the detection of particles with energies 
in excess of 10 20 eV. Acceleration processes are 
strongly limited by energy losses and finite size 
of the known acceleration regions and only some 
types of sources are barely able to energize pro- 
tons up to the observed energies pq |. Compo- 
sition and anisotropy studies will be the keys to 
solve the mystery, but at the cost of increasing 
the statistics by at least a factor of 10 compared 
with current experiments. Two experiments are 
being planned for the next decade or so, and 
will provide the characteristics necessary to do 
cosmic ray astronomy: the Auger project f37|| is 
currently in the construction stage in Argentina, 
while the EUSO project Q is scheduled for op- 
eration starting in 2008. Each one of these enter- 
prises implies an improvement by a factor ~ 10 
compared with the previous one, which means a 
predicted 500 events per year above 10 20 eV for 



EUSO, if the AG AS A spectrum is taken as a tem- 
plate. 

Besides being the tools for ultra high energy 
cosmic ray astronomy, these experiments repre- 
sent a unique tool to study possible New Physics 
at extremely high energies. The case of neutrino 
oscillations provides an example of the first hint 
of the existence of Physics beyond the Standard 
Model of Particle interactions, derived in an As- 
trophysics context. It is foreseeable that the ball 
of particle physics, after a few decades, could go 
back to the field of cosmic rays, where the first 
steps in that direction were moved in the '30s. 
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